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FOREWORD 


The work reported herein, covering the period from 1 November 1969 through 
31 December 1972, «u'carried out by the Infrared and Optics Division of the 
Willow Run Latx>ratones, then a unit of The University of Michigan's Institute of 
Science and Technology. (On 1 January 1973. the Willow Run Laboratories sepa¬ 
rated from the University and became independent as the Environmental Research 
Institute of Michigan iERIM), P.O. Box 618, Ann Arbor, Ml 48107.) The Air Force 
Avionics Laboratory (AFAL) of the Air Force Systems Command, Wright-Patterson 
Air Force Base, Ohio, commissioned this work under Contract F33615-70-C-1123, 

Project 6239, Task 10. Mr. Bruno Wernicke AFALv'RSP is Technical Monitor for 
the Air Force. 

Oir interest in the development of a bidirectional reflectance model is an out¬ 
growth of on-going research In target signatures which has thus far produced eleven 
Data Compilations listing reflectance characteristics for materials of interest to 
the Target Signature Analysis Center (TSAC) and its patrons. 

Some work toward extending the mooel to account for a non-Lambertian, noo- 
specular reflectance component assumed to result from scattering within the target 
material was performed under Contract DAAD05 -7J-C-0246 with the Army Ballistic 
Research Laboratories (BRL). Under the same BRL contract, the model was also 
coded in Fortran IV and an extensive measurement program implemented and used 
in a validation of the complete model. Because the present unified model nas been 
realized as a result of work- performed at WRL/ ERIK under the present Air Force 
Avionics Laboratory contract as well as from work done under the BRL contract, 
the BRL-sponsored portion of the model is also described in this report. 

This research effort is continuing at ERIlf's Willow Run faculties under the 
direction of Dr. Robert Maxwell as Principal Investigator, with guidance provided 
by Mr. R- R- Legault, Director of the Infrared and Optica Division. The ERIW num¬ 
ber for this report Is 196400-1 -T* 

ThU report was submitted by the authors on 11 July 1973. 

This technical report has been reviewed and is approved. 

Albsrt W. Berg, Chief 

Reconnaissance Sensor Develo pme nt Branch 

Reconnaissance Division 

Air Fores Avionics Laboratory 





ABSTRACT 


This report describes a method for using bidirectional reflectance information 
previously reported in the Eleventh Supplement to the Target Signature Analysis 
Center; Data Compilation [ 1, 2] and further validates the bidirectional reflectance 
model originated and extended under recent contracts. It includes bidirectional re¬ 
flectance model parameters for a variety of paints. Parameters -were extracted 
from measurement data reported in the Eleventh Supplement . Reduced reflectance 
data are also provided: these data may be used with the computer model or option¬ 
ally. in an interpolation procedure for estimating reflectances without the aid of a 
computer. 

The computer model makes it possible to calc-tale !'directional reflectance data 
from a very small amount of measured data. Accuracy demonstrated in the Model 
Validation section indicates that the model is very effective, although improvement 
can still be obtained at large receiver zenith angles. The interpolation procedure 
also shows excellent agreement with measurement. 
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BIDIRECTIONAL REFLECTANCE MODEL 
VALIDATION AND UTILIZATION 

1 

INTRODUCTION 

A model for predicting the radiance at a remote sensor must include the spatial, spectral, 
and polarization characteristics of the bidirectional reflectance and directional emittance with 
respect to target and background surfaces. In prua:iple, the directional reflectance and di¬ 
rectional emittance properties of materials must be known for all source and receiver angles, 
polarizations, and wavelengths. A Lambertian assumption may be valid for some types of 
backgrounds, but for most man-made targets is scarcely adequate. Measurement of all spatial, 
polarization, and spectral characteristics of the bidirectional reflectance and directional emit¬ 
tance for a large number of material samples is impractical. Even if such measurement were 
performed, the data could not all be stored efficiently enough to make it accessible for digital 
computations. Clearly, an empirical model is required to approximate the ••■.directional reflec¬ 
tance and directional emittance properties from a limited number of measurements. 

The bidirectional reflectance model developed by the Environmental Kcsearch Institute of 
Michigan (ERIM) is described in this report. The model accounts for effects that produce both 
specular and diffuse components. In pertlculai, a surface model relates bidirectional reflec¬ 
tance for all source-receiver angles and polarizations to fixed bistatic measurements and a 
Brewster angle measurement. In addition, the model enables calculation of either a Lamber.ian 
diffuse component or a non-Lambertian diffuse component. The Utter component accounts for 
anguUr and depolarization properties arising from internal scattering effects. Our extension 
of the bidirectional reflectance model has considerably improved the fit between model predic¬ 
tions and measured data, as will be shown In Section 0. 

( 

In addition, we calculated fix tit' bistatic data from reflectance data for 20 material# included 
in the Data Compilation of the Eleventn Supplement [2]. Results of these calculations are 
graphed and tabuUted in Appendix I. 

A useful method lor deriving reflectance data when no computer Is avalUble is also pre¬ 
sented. For this purpose, simulated data representing typical sets of parameters are provided 
as well as a method whereby typical fixed bistatic data may be used to bracket measured zero 
bistatic data for a particular material so that the bidirectional reflectance for that material 
may be estimated by an Interpolation method as described In Section 8. 
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As l*. now stands,the model permits generation of an enormous amount of bidlrect.onal re¬ 
flectance data from a very snail amount of measured data. The accuracy shown in Section 6 
on Model Validation indicates '.hat the model is very effective, although it can still be improved, 
particularly at large receiver zenith angles. With the abti'ty to account for elliptical (partic¬ 
ularly circular) polarization now built In, the model is available for use with circularly pc - 
larized sources, if these sources prove useful in the future. 

In this report, we compare measured data with results computed from both the initial 
model and from the extended model,and then evaluate the relative performance of he two, W’e 
establish a domain of validity for each, based on material properties. Since the modeling is 
empirical, only a limited amount of measured data are required as input parameters. In thie 
ease, the parameters are the fixed bistatic data. Therefore, data from the Dnta Compilation 
[1, 2] are used to derive fixed bistatic curves contained in Appendix I. 

All modeling described in this report was performed with respect to oot wavelength, A - 1.06 

M ni. 


a 







BIDIRECTIONAL REFLECTANCE 


One pnysiral property which can be measured directly from a sample of material is bi 
airectional reflectance. The physical definition is 


p'H .0.; - .0 1 
i 1 r r 


:L 1“ .) 1 
r r r 

'£!■■, 3 j 

i i l 


(la) 


where C/E.l-)., ? ) is the incremental irradiance (power per unit area) impingent on the surface 

of a material frornthedirectionf-,. i ). and L (■* , b ) is the resulting increment of radiance 

i i r r r 

(power per unit projected area per unit soiid angle! scattered from that surface in the direc¬ 
tion (- . * r ). Figure 1 illustrates the situation. The bistatic ancle. £!, is that angle between 
the vectors which point to the source and the receiver respectively. 


Equation (la)canbe rewritten interms of directly-act -ssible experimental parameters as 
SP 
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where c p is the power, in watts, incident from the direction <•> . * ) or. the small area CA, and 

i ti 

is the resulting power scattered Into the small solid angle Ofl r In the direction (o^). 

When polarization dependence is to be shown, subscrip are appended to thep' term. Thus 

when we wrltep’ . the leading subscript, a . describes the source polarization while the trail¬ 
er 1 

lag subscript, a^. describes the receiver polarization. The source polarization, always re¬ 
ferred to the plane of incidence, describes the polarization state of the electric field vector. 

The appended subscript symbols I and I indicate whether the source electric vector polariza¬ 
tion Is para;) i o or perpendicular to the incidence plane. The reflected electric field polariza¬ 
tion state tt xpictfied by the same symbols, but here the reference plane is that ref'ertance 
plane defli /. try the sample normal and the direction to the receiver. (For example, o'. re- 

A ■ 

present«. ,, Recta nee messsred when source pofc-rtzatton is perpendicular to the incidence 
plans as./ receiver polarisation la parallel to the reflectance plane.) Notice that when either 
the source or the receiver, or both, are scanned In angle over the sample, the tncldenre and 
reflectance planes change orientation with relation both to the sample and to each other. 

Bidirectional reflectance depen ds on the physical properties of the material aa well as on 
the geometric state of its surface. Different aortic-' states result In different reflectances. 
Hence, a complete collection of bidirectional reflectance data (or any single mats.-lal would 
require measurements of a large number of samples of the matertaLewch with a different sur¬ 
face state. Each sample would have to be measured with several soerre-reeeiver polarization 
comMnations. Con sequent ly a very large number of source and receiver poeitione would be 
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required (or each set of polarization states. Finally, the enure procedure would have to be 
repeated at many different w-aveleneths. The purpose of modeling is to predict reflectance 
data from only a limited number of measurements and hence eliminate the need for an other¬ 
wise unwieldy measurement program. 
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BACKGROUND INFORMATION 


The bidirectional reflectance model Is based on observations of polarized bidirectional 
data from rough, painted surfaces which exhibit a Brewster angle (Fresnel-like behavior in 
relation to the specular geometry). The degree of depolarization appears slight, and oc that 
basis (or the initial modeling work in this program single specular reflection from the rough 
front surface was assumed to be the dominant reflection mechanism. Multiple front-surface 
reflections and internal scatterings were observed to be smaller and were initially incorpo¬ 
rated into a Lambertian "volume" model to account for the diffuse component. 

The assumption that the diffuse component Is Lambertian, however, makes it difficult to 

account for certain anomalies that occur when measured data are compared with the model's 

output. For example. Fig. 2 is a bidirectional reflectance curve showing the reflectances at 

the receiver as the receiver scans over zenith angles from 0° to 90° in the ♦ =0° and * = 180° 

r r 

hall planes. The source remains fixed at A * 40° and e » 18CP. The upper curve shows re¬ 
flectances when source and receiver are both linearly polarized at the same polarization angle 
with respect to the target-Incidence and target-receiver planes. (In this case, both are per¬ 
pendicular -polarized.) The lower curve shows reflectances when source and receiver are 
cross-polarized with respect to one another. (Source is perpendicular-polarized; receiver is 
parallel -polarized.) Note the marked angular d ep e ndence la the lower curve. If the nonspecular 
c o mp onent were truly Lambertian, no such angular dependence would be present. 

Also, although -adlation sources in this work are all linearly polarized, future work may 
well Involve more general cases. Therefore, the model should account for the most general 
type of polarization—namely. eUlpticaL 

For the above reasons, and In order to obtain a closer overall correspondence between 
model prediction sad means red data, the model has been extended to account for the following: 

(I) possible ncm-Lamb*man angular dependence of depolarised component 

(X) shadowing sad obnea ratios produced by the roughness of the surface 

(3) elliptical polarization 

The model—a pheaomssological one in that Its sen required a limited dumber of measure- 
meats —id described is the nest two sect i o n s. Section < vc lodes a discuss ion of specular re- 
flsctaace from thd surface, effects canned by shadowing anu obacsration resulting from surface 
roughness, and polarisation effects. Section S describes the volume model. 
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SURFACE MODEL 


!jj this section, we review the surface model and also discuss the interference effects that 
necessitated model modification. 


4.1. AVAILABLE AREA 

If the rough surface is considered to be made up of small sequins having a distribution of 
orientations, there will be some specular reflectance at any receiver angle and the extent of 
that reflectance will be determined in part by the amount of surface oriented for specular re¬ 
flection at that receiver angle. (The area available for such reflection will also depend on how 
some sequins shadow or obscure others.) Since measurements do, in fact, show a reflectance 
distribution over the hemisphere, we assume that the above description is valid and that there is. 
Indeed, a distribution of surface areas which have normals pointing in different directions. There¬ 
fore, to establish the distribution of available surface area, we define a density function Z(8, 4) 
which describes the relative density of local surface normals (per steradian) pointing in the 
direction ( 9, 4). 

The effect of the distribution of surface normals is measured by a zero btstatic measure¬ 
ment in which 9 - « r and * * r - (Note that we really use a fixed bistatic scan with a small 
bistatic angle, A true zero bistatic scan wrold be very difficult to obtain since source and 
receiver obviously cannot occupy the same position.) 


4.2. FRESNEL COEFFICIENTS 

Fresnel reflectance coefficients describe the reflectance and polarization of specularly 
reflected radiation as functions at source and detector positions and of the complex index of 
refraction. However, since we are trying to find reflectance as a function of source and de¬ 
tector positions only, we must know—or he able to determine —the index of refraction. (As 
(Harassed later in this section, if can determine the index by measuring the Brewster argle.) 
Since, la the surface model, we consider only single, local specular reflections, the Fresnel 
equations automatically account for polarization. 


IT the receiver subtends tbs solid angle 90^ from tbs s am pl e (ess rig. 1) tbs solid angle 
6Q* la which local surface normals must Ue to permit collection of tbs local specularly re¬ 
flected radiation by tbs receiver la gives by: 






This solid angle la c entere d about the direction 



( 2 ) 











Moreover, k is taken to be very small* so that n can be determined experimentally by mea¬ 
suring the Brewster angle, 9_, and then using n = tan 9_ to solve for n. 

fi D 


4.3. SHADOWING AND OBSCURATION 


Referring to Eq. (9), we can derive a zero btstatic curve, p ( 9 . ,0 A ; 0 A ,* A \, from a pt0,d., 

\ n n n a/ it 

9 4 } curve with 9 , ♦ fixed and 9 variable by inverting the equation so that 


P '{9 d , 9 ♦ 


,p‘(9.*.9 d ) cos 9. cos 9 


ft a' a a/" 


cos 9 


after doing this for a variety of 9's. we found that the curves obtained differed systematically 
from those obtained from a fixed bistatic measurement. Apparently, because of surface rough¬ 
ness, some sequins shadow or obscure others; this reduces reflectance everywhere except at 
a purely back-scattered position. The model must therefore be modified to correct for such 
Interference. 

Torrance and Sparrow [4] have developed an analytical function that helps correct the sit¬ 
uation; however, we have constructed our own function using empirical considerations only. 

Oir function results in better agree m ent between measured and derived fixed bistatic curves 
than does the analytical function at Torrance and Sparrow. The empirical function (SO) la de¬ 
fter i as: 


»♦*•** ( 

A ' 


V 1 


where O sad r are parameters, and is s factor calculated from the geometry, which ad fusts 
the fall-off rate of the shadowing and obecaration function tat the forward-scattered direction. 


•For the calculations in hits study, results of past measurement programs [1] were 
seed to establish the refractive Indices. In those programs. It was determined that the 
magnitude of the total index of refraction was close to 1.03; that the imaginary part of the in¬ 
dex of refraction could be seglected, compared to the real part; and that the index of refraction, 
for the wavelength* of Incident radiation under constderattoa (1 to 4 pm), did not vary apprecia¬ 
bly. 
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VOLUME MODELS 

The following discussion outlines the reasoning behind the extended portion of the bidirec¬ 
tional reflectance model. (The extended portion Is referred to as the "volume" model.) 

Different materials with varying degrees of surface roughness and different optical prop¬ 
erties show differences in nonspecuiar reflectance behavior. These differences show up in the 
extent to which the nonspecuiar reflectance is dependent upon angular position of the receiver. 

To make provision for materials that do exhibit such angular dependence and for those that 
do not, two volume models are used. The following discussion describes, first, a Lambertian 
volume model which has no angular dependence, and then a non-Lambertian volume model in 
which angular dependence Is Important. 

5.1. LAMBERTIAN 

In addition to Fresnel reflection from a surface, other effects sum as might take place be¬ 
neath the surface can produce a nonspecuiar reflectance component everywhere in the hemi¬ 
sphere. if the surface roughness as well as the absorption propertir s of the surface are right, 
this volume reflectance may be completely diffuse and uniform over the hemisphere. Moreover, 
the reflected radiation will be totally depolarized, regardless of the polarization of the source. 
Than, If the receiver Is polarized In the orthogonal direction to the source polarization, an ln- 
ptane measurement will represent the volume component only. However, only half the volume 
c o mpon e n t Is actually represented, since there should be an equal diffuse contribution polarized 
In the same direction as the source. 

The Lambertian volume component is one of the input parameters for the model when a 
target material with Lambertian reflectance properties Is considered. A method whereby values 
for this parameter may be extracted Is described in Section 6. 

5.2. NON-LAMBERTIAN 

Ob the basis of the Lambertian diffuse model described above, no angular dependence 
would be expected for the diffuse component. However, for some materials, actual measure¬ 
ments show that there Is an angular dependence. To provide for the angular dependence of the 
diffuse co mp on en t, the model has been extended by Including scattering that takes place beneath 
the surface. 

Assuming an exponential scattering function as the radiation first enters and then leaves 
the surface, and making reference to Pig. 3, we construct an expreeaion for the volume scat¬ 
tering component of the bidirectional reflectance as follows; 


















E iO = Irradiance at surface of area A^, where A^ Is the area of cross section of the colli¬ 
mated beam and is normal to the beam 
Eg = irradiance on surface element of area A 

= irradiance on surface of slab or area A at distance a beneath surface 

L r (0) = radiance scattered from primary beam through 20 In direction of receiver 
0 = half of angle between target-to-source and target -to-reeetver vectors 
a * total scattering cross section (ignoring absorption) 
d(fl) = differential scattering cross section with respect to 0, l.e., fa(0)dQ~ Jdo/dftdft =<r 
Q * solid angle subtended at target by receiver, assuming a point target 
^ » angle of Incident beam relative to fixed z axis 
* angle of reflected beam relative to fixed z axis 


The objective of the following calculation is to determine that portion of the primary beam 
scattered from distance z beneath the surface through an angle 20 toward a receiver which sub¬ 
tends solid angle ft. 


First, the bidirectional reflectance defined In Eg. (1) is now o' = Lj/Eq with respect to the 
slab (see Fig. 3). To determine E Q : 

A^ * A cos 8^ 

P/A, 


p 

E„ => t 


P cos 


0 A cos 6. 


E i0 cos 8 { 


i A I 

where P Is the power at surface of ares. A. 

The irradiance incident on the slab at distance z beneath its surface is: 


. -oz/cos 6. 

Ej » EgS -01 = Eg® 1 


where f » x/ cos 8^ Hence 


-or/cos 8. 


*1 * E 10 c0 * *l e 


1 


and 


-ds/cos 8, 


* - E io°* 


(14) 

(15) 


(16) 


( 17 ) 


(It) 


where dE. Lb the amount by which Irradiance decreases In going from distance s to distance 
1 -dt/eoe 8 t -d s/cos 6, 


s ♦ dz beneath surface. Note also that e * and • 1 represent the scattering 


loss from the beam on the way In and on the way out of the material, respectively. To deter- 









In Ignoring the finite thickness of the layer of material, we have also Ignored the possible 
specular reflectance of the bottom surface. To account for the possibility of specular reflec - 
tion from the bottom layer, it may be useful to provide a parameter function peaked near 
9 s * 0. Therefore, we include all 3 dependence in a function 10), and all 9 A dependence In a 

a / \ n 
function g|0 A j, and write 


P* 


2-5_i_ 

COS + COS 0 r 


(28) 


where f(3) and g($ A [provide freedom fcr empirical adjustment. The constant, p^ represents 

/ \ 

the value of p' when = 0 and 10) = g(d A j = 1. 
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MODEL VALIDATION 

Use of the bidirectional reflectance model requires a limited amount at measured data 
namely the zero butane measurement) from which complete sets at reflectances can be- cal¬ 
culated The results of these model-calculated bidirectional reflectances can then be com¬ 
pared to corresponding resuits of actual measurements. This was the procedure we followed 
to validate the model. 

Model calculations and measured data were compared sn terms of p i the reflectance), a or 
„ ^ (the ancle of polarization for ;r,«* oeam alter reflection from the target), and P (the percent¬ 
age oi polarization of the reflected beam). 

Measured dais for materials at different properties (color and roughness) were used to 
demonstrate the model's performance. The materials are designated as A030I8-00I and 
A02018-002 and are not included in the previous compilar.on j 1.2|. Material AQ 2018-001 Is a 
green paint and material AO2018-002 is a tan paint. These materials were supplied by the 
Army Ballistic Research Laboratories for the purpose of dsveloping the non-Lambertian dif¬ 
fuse component cf the model. They were used for validation because an unusually extensive 
set of measurements could be made on them. Model parameters are listed in Table L (See 
Section 7 for definitions of model parameters.) The overall discussion cf the model fitting is 
divided into four parts: 

(1) p for A02018-001 

(2) p* for A02018-00J 

(3) polarization angle (a or lot A02018-00! 

(4) percent polarization (P) for AA201I-001 and A02f It-003 

tn what follows, the or natation cf the source polarizer in the measurements of materials 
A02O18-OOI and A02Old-OO3 was sot actually perpendicular, parallel, sor at 43° to the plane of 
incidence but Instead was offset by 3° In sack ease. Specifically, the appropriate correspon¬ 
dences, s hown in Table D, should be recog n ized. These shifts were taken Into account when 
the validation calculations sere mode on the computer; however, we cost le ws to refer to ‘'per¬ 
pendicular,” "parallel." and 4$°~ 

II. REFLECTANCE FOR SAMPLE MATERIAL A03011-001 

Material ACrJOll-OOl is a green painted surface. The aero Mstsilc measuremen t with 3° 
polarisation angle (I#., alasost perpen di cular polarisation) ta shown tn Fig. 4. (The aero bt- 
stafic data w«h parallel -polarised source, although not shown, have Identical characteristics.) 
The aero b<static plot ta sh ar pl y pee k e d cf 0°, falling off rapidly to a constant rales at about 


IT 








TABLE I. MODEL PARAMETERS FOR SAMPLE PADfTS 


Material 


Parameter 

A02018-001 

A02018-002 

n 

1.65 

1.65 

k 

0 

0 

p xl 

_* 

0.044 

p x2 


0.044 

P 

0.007 

0.05" 

V 



T 

15 

15 

a 

40 

40 

( 03 ) 

1 

1 


1 

l 

\ <»/ 



a ’& • B 9 ■ 0 

K a’ A' A f A , A 



\ n n n n. a 



X 

1.06 fia 

1.06 pn» 


•Thia material l> ran with ’!>• noa-Lambertian rota me model: therefore 
▼ aloes are not aKeuary. 


''Material 2018-002 was rw» with the Lambertian rotume model: there- 
(ore Py ahoold not be need. 


TABLE O. TRDI SOURCE POLARIZATCft ANGLES 


rtlCtiftf 


Woml—1 Angle* 


Act met) 
Plane 

8 

1P°1 

1(60°) 

-ft® 

-45® 

0°-lS0°' 
S0 # -210° 
•0°-240°, 

\ 

S° 

-«® 


-46® 

•0 # -270® 


S* 

*66® 

♦50® 


rtaed Metatic 

5° 

•65° 

»W° 
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20°. In all receiver polarizations, p' shows an angular dependence clearly departing from 

! \ 2 

Lambertian behavior. The table, at values for p ^ „ )cos e . used in the model 

' h n n n' n 

wan obtained from this measurement by reading ctf p‘ and p. at each angle and. then calcu- 
2 *' 

Laling ip' - p' ) cos e where 9 is the angle that the normal to the reflecting facet make# 

i > • 1 • A A 

* 1 a 3 

with the fixed z-axis. In zero bistatic scans, 9^ = 9^ = 8^ (see Fig. 1). (Physically the source 
and receiver were separated by 1.8°. Thus, both were 0.9° from the true 9 . In the calcula- 

A 

B 

tions, the axis was translated to bring the x-axts into correspondence with t) =0.) The sub- 

n 

traction, p| - p' ( eliminates the diffuse contribution which would distort the value for 
a'U ,6 i which is what must be measured (recall Eq. 9). 

A A A *. > 

\ a n n n' 

o 

In Figs. 5, 7, 9 and 11,* plots of measured data are shown for » = 40, ♦. = 180 and where 
is scanned in azimuth planes represented by = 0°, ISO 0 ; 90°, 270°; 20°, 210°; 60°, 240°. 
Each measurement plot is followed by plots ct data generated, respectively, by the Lamber¬ 
tian model with no shadowing and obscuration factor, by the non-Lambertian model with no 
shadowing and obscuration factor, and by the non-Lambertian model with the shadowing and 
obscuration factor. For example. Fig. 6 shows the calculated p data for ^ = 40° and <* r as 
scanned in the 0° and 180° azirmxh planes for the above variations of the model. The simu¬ 
lated source is taken to have s 'pe r pendicular'’ polarization angle. In these in-plane scans 
(♦ r * 0°, 180°), the main peak is ui the 0° azimuth plane which is the forward direction for the 
source angle at ♦. * ♦ 180 °. Nate the rise (in the plat at measured data) at large zenith angles 
for the cross-polarized component. This is a characteristic which suggests the coed for the 
non - Lambert tan Volume model. 

Surface Plus Lambertian Volume Model with No Shadowing and Obscuration Correction. 


Figure • plots (in solid lines) the model calculation using the surface model plus the Lamber¬ 
tian volume modal with no correction for shadowing sad obscuration. Tbs following charac¬ 
teristics should be noted; 

(1) In the 4 • 0 (forward scattering) azimuth plant, the modal fits the measured data 

very well b etw een * f • 0 sad * f • S0° far marcher! polarization of source and receiver. 
At 0 f - SO®, the calculated curve suddenly diverges. This Is thought to be the result 
at the failure to acc oew t for Aadowtng and obscuration as di s cuss ed earlier. At # f 
• 0° anti on into the beekscaltered (# r * ISO 9 ) direction, the calculated values lie 
above the measured values and this, too, la believed to be the result at the tack at a 
shadowing and obscuration correction. 


(2) in the cross-poUrtulion component (il), the model predicts a flat response except 
for a alight hemp under the specular peak. Tbs measured data, koeeeer, show a 
Osar angular d epends ne e ee # f . 


ate: Oh ail reprints at original computer plots, the symbols $ sad d are represented 
by sod respectively. r r 
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With the exception of these two characteristics, however, the surface plus Lambertian volume 
model with no shadowing and obscuration correction fits the measured data fairly well. 

Non-Lambertian Volume Model with No Shadowing and Obscuration Correction. The dotted 


line* in Fig. 6 show a mode' plot using the tame parameters, except that the non-Lambertian 
volume model is now used. Keep in mind that one may use the non-Lambertian volume scat¬ 
tering as a model by itself or in conjunction with a specular component. The latter Is used 
here. In the like-polarized component, nothing has changed from the previous case. However, 
the cross-polarized component now fits the measured data much more closely. It rises steadily 
at large angles, both in the beck-scattered and forward-scattered directions—a result of 
L (cos d + cos 6 f ) dependence shown in Eq. (28) for the volume model. However, the response 
for the like-polarized component does not drop sharply enough at either side of the peak, and at 
high angles in the forward-scattered direction, the awkward divergence still appears at 60°. 

Thus, the non-Lambertian volume rvrtel i mpr oves the cross-polarized fit (with respect to 
material AQ2018-001) over that of the Lambertian volume model and, apart from anomalies at 
high angles and near 0°, provides a reasonable fit to the measurements. 

Non-Lamberttan Volume Model with Shadowing and Obscuration Correction . The dashed- 
line curves in Fig. 6 show results with the shadowing and obscuration correction applied to the 
nos-Lambertian volume model calculation. The cross-polarized component is unaffected. The 
net effect an the match-polarized component is to reduce the reflectance everywhere except at 
the specular peak and at the direct beckacattertr* peek (i.e., at 0 * 0). In particular, it lowers 
8 - forward-scatter contributions beyond 50°, bringing the model closer to measured data in 
this region. Overall, the fit obtained using the volume model with a shadowing end obscuration 
correction agrees closely with measure menu. 

The forego in g discussion applies to "in plane" receiver scans—those in the * 0 and 
■ 100° azimuth planes. The azimuth plane perpendicular to the 0°, 180° plane is the 90°, 

270° plane and is ref e r re d to as "out-of-plane". The plane we are in or out of Is the plane of 
incident beam and target normal, or the target incidence plane. (See Fig. 1.) 

Is Fig. 7 we have the plot of measured data for the out-of-plane situation with perpendicular* 
polarized source agala. la this case, ho wever, the Incidence plane Is per pe ndicul a r to the re¬ 
flection plane. At ■ 0, therefore, p|jte plane la the same ae p| j otd of plane, for this 
reason, the reflectances of match-polarised and croee-polartzed comp o nen ts seem to exchange 
behaviors In the out-of-plane configuration, as Is verified by the plotted measurements as well 
ss by the model calculations. Figure • presents {dots of a Lambertian model without the shad¬ 
owing and obscuration factor, a surface plus non-Lambertian volume model without tbs shad¬ 
owing and obscuration factor, and the surface plus aon-LambertUa volume model with the 
abed owing tad obscuration factor. Ae before, It tr ap pa r e nt that the use of the noa-Lanfeertlsa 
votes* model plus the shadowing and obs c ur at ion factor Improves a gree m eat between modal 
and mea surem en ts so Out, apart froaa a possible overall scale factor, the agreement Is within 
measurement fluctuation. m 
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For additional validation, plots are shown for the 30°, 210° azimuth planes (Figs. 9 and 10) 
and for ‘.he 60°, 240° azimuth planes (Figs. 11 and 12). The characteristics of calculated and 
measured curves, apart from a scale factor, are in excellent agreement. Figures 13 through 
20 represent similar comparisons for the case when the source polarizer is set for -45° (In ■ 
the Q v , 180° 'tzimuth plane) and set parallel (In the 30°, 210°; 60°, 240°; and 90°, 270° planes). 
Measure:, plots are presented with the calculated plot to represent the surface piss non- 
Lambertian volume model and to include the shadowing and obscuration factor. 

6.2. REFLECTANCE FOR SAMPLE MATERIAL A02018-002 

Material A02O18-OO2 consists at a tan painted surface. 

Based on the zero bistatic scan, Fig. 21, material A02018-002 appears to fct- somewhat 
brighter than material A02018-001. Whereas the non-Lambertian volume model was clearly 
the best choice for material A02018-001, it is not in the case at A02O13-OO2. In this latter 
case, the best choice Is the Lambertian model. . . 

The lack at angular dependence In the reflectance of the cross-polarized component could 
have a number at explanations. Multiple scattering increases for rougher surfaces. Since such 
scattering may nek be angular dependent, it could become a large enough factor to swamp the 
angnlar dependence which is otherwise present. Moreover, the difference in color b etwee n the 
green and tan certainly alters the absorption sad, consequently, can alter the angular dependence 
as well. 

In any case, the appropriate model to use can be determined by looking at the cross- 
polarized component at the fixed bistatic scan, ft a clear angular dependence is present, the 
non-Lambertian model should be used. But if there Is little or no apparent angular dependence, 
as with material A02018-002, then the Lambertian model Is more appropriate. 

b Figs. 22 through 29, plots are p rovided for different azimuth planes, beginning with the 
plot for measured data, followed Immediately by the corresponding plot from model calcula¬ 
tions. In this grcxip at illustrations. Fig*. 22 through 25 re p re sen t pe r pend i cular source po¬ 
larization, while Figs. 28 through 29 repre s ent a scarce parallel polarization for the 0°, 180° 
axtsndh plane and for the 90°, 270° azimuth plane. 

In all cases the fit appears to be excellent, except for occasional anomalies at targe sal- 
math angles. Further modification at the shadowing and obscuration factor should decrease 
these present anomalies. - 

8.3. POLARIZATION ANOLT WJTOR 8 AMP LX MATERTAL .*02018-001 

The reflectances of the perpendicular tad parallel cowpfcs'.-Ma of a linearly polarized 
beam vary as functions of the soarce-recetver angles und the Index of refraction of the target 
mafertaL (See, for example, the Fresnel equations, Paf. 3.) Based an observations, the index 
of refraction varies little over a wide range at paint surfaces. For the parties tar materials 
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FIGURE IX CALCULATED p' FOR A0301B-001 USING NON-LAMBERT IAN VOLUME MODEL 
WITH SHADOWING AND OBSCURATION FACTOR. 9, « 40°, p, « 180°, * - 60°, 240®. 
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covered in this report, the variation is considered to be zero. Therefore, for these surfaces, 
polarization angle is essentially a function cr.iv r{ source-receiver positions. 

As additional validati.cn for the model, predicted polarization angles are compared with 
polarization angles extracted from the measured data. Figures 30 through 33 show plots 
obtained for the O' 3 , loO°: 30', 270°: 30 3 , 210 and 60°, 240° azimuth planes. Measured data 
represent material AC 2018-001. In all cases, agreement between measurements and calcu¬ 
lations is excellent, with the average disoarttv not more than 10^-. In particular, the dramatic 

Q O OO 

agreement between measurements and model m the 30 , 210 and 60 , 240 azimuth planes 
constitutes powerful verification at the mcoel and affirms its usefulness in arbitrary source- 
receiver •xieittcws. 

6.4. PEPCENT POLARIZATION POP SAMPLE MATERIALS A02018-00* AM) A02018-002 
Percent polarization ■ P) validate* the ratio of surface-to-volume contributions to reflec¬ 
tance. Percent polarization depends on both polarized reflectance and angle of polarization, 
both validated in earlier sections d this report. In this section, we compare mode) predic¬ 
tions with percent polarization values extracted from measured oata. 

Figures 34 and 35 illustrate degree of polarization for scans at material AC2018-002, for 
p erpe n dicular and parallel sources, respectively. The vaiidifv a t the model Is supported by 
the clone correlation between the behavior e< values extracted from measured data and those 
calculated with the model. 

Additional confirmation of the model Is -.rovtded In Figs. 36 through 38 where percent 
poianzjtian plots are given for material AMO 18-001 la the 0°, 180° and 90°, 270° azimuth 

angle planes. 





















O: Measured Source, *5° 
▼: Measured Source, -40® 
a: Measured Source, -85® 
x: Calculated Source, t-3 . 
©: Calculated Source, -40'"' 
A: Calculated Source, -85' 
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FIGURE 34. PERCENT POLARIZATION 
VARIATION FOR AQ2018-002 AS FUNC¬ 
TION CF SOURCE-RECEIVER POSITION. 
3 = 40°, 6. * 180°, * = 0°, 180°; perpen- 
1 dlcular source. 



FIGURE 35. PERCENT POLARIZATION 
VARIATION FOR A02D18-002 AS FUNC¬ 
TION OF SOURCE -RECEIVER POSITION. 
6 =. 40°, * 180<>, « r - 0°, 180°; parallel 

source. 
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FIGURE 36. PERCENT POLARIZATION 
VARIATION FOR A02018-001 AS FUNC¬ 
TION OF SOURCE-RECEIVER POSITION. 

9 • 40°, * t . 160°, 6 r » 0°, 180°; perpun- 
dMuUr source. 


- MX’ <■*')■■ 

























MODEL PARAMETERS 

This section briefly describes the model parameters that can be used in the bidirectional 
reflectance program and explains how their values a-e derived. The choice of parameters for 
use in the program depends to some extent on the mode of the model being run. Basically, the 
model is run in three different modes: 

(t) Surface and Lambertian volume components 

(2) Non-Lambertian volume component (no surface contribution included) 

(3) Surface and non-Lambertlan volume components 

Therefore, we have grouped the model parameters as follows: 

(1) Polarization parameters 
{2) So rface model parameters 

(3) Lambertian volume model parameters 

(4) Non-Lambertian volume model parameters 

(5) Parameters used to generate typical data for comparison purposes (see Sec. 8) 


7.1. SOURCE POLARIZATION PARAMETERS 

The present model has been designed to account for polarization dependence in both sur¬ 
face and volume components. 

In the surface component, polarization is accounted for automatically In the Fresnel reflec¬ 
tance coefficients. In the most general case, such polarization can be elliptical and can be de¬ 
composed Into linear and circular components. To date, only a linearly polarized source and 
receiver have been used In measurements. However, for some applications, circularly polar¬ 
ized sources or receivers may be of interest. Therefore, In the model, we have provided pro¬ 
gram subroutines which take into consideration the eiliptlcity and handedness (t.e., direction 
of rotation In an elUptlcally polarized source) of both Incident and reflected beam. 

For volume components in both Lambertian and non-Lambertlan cases, it Is assumed that 
reflectance will be depolarized to some extent. In both cases, In fact, we assume total depolar¬ 
ization. Therefore, although a depolarization factor has been included In the non-Lambertian 
volume model for future flexibility, we assume DP(0) * 1. 

TS source polarization may most generally be defined ss partially polarized wkk the po¬ 
larized component elUptlcally polarized. The stats of polarization of the source wtU be defined 
by Its degree of polarization, P, and parameters A, B, if, and H to define the elliptical polariza¬ 
tion of the polarized component Hare, A and B ars the intensities along the semi-major and 
semi-minor axes, respectively. The angle ^ Is the angle between the semi-major axis of the 
elUpse and the direction normal to the plane of incidence, measured looking Into the source *• 
beam; if Is equivalent to a except that 0° s i> s 180° and -90° sa* 90°. The handedness H 
• *L ‘ ' 














The Stokes vectors pxovide a convenient formalism for defining the polarization state of 
the reflected radi*.i.ce. (Reference [5] proviaes a general discussion of Stokes vectors in this 
context-1 


V‘u 

I cos 2x cos 24\ 

P 

I cos 2x sin 24 
P 

1^ sin 2x 

where ! and are the polarized and unpola'rized components, respectively, in the reflected 
radiance. The degree of polarization in the reflected radiance is P = I+ 1^). Angies x 
and 4 define the polarization state of the reflected radiance; 4 is the angie between the semi¬ 
major axis of the ellipse and the direction normal to the plane of reflection; tan x = aVB/A 
where A and B are the intensities along the semi-major and semi-minor axes of the polariza¬ 
tion ellipse and tan X ^ q for handed elliptically polarized radiation- 

The RHOPR1ME program produces the Stokes vector S for unit irradiance in the input beam; 
the area may also be defined to be unity and then S represents a reflectance Stokes vector. The 
program also produces the components of the reflected radiance transmitted with a receiver 
polarization analyzer oriented parallel or perpendicular to the reflectance plane for computing 


7.2. SURFACE MODEL PARAMETERS 

One of the quantities in Eq. (9) from whichp'(0„*.; 0,4 )ls determined isp'le.^ ; 8 ,6 \ 
* , \ a h n n / 


cos* f,. A» previously discussed, p’f& A ,4 A : 

2 ' 


A" A’ ' ^a)“ 


obtained from zero blstatic data. Values 


.an a 

for p'/lf ,♦ ; 6 ,4 \coe 8 must be calculated (preferably for Increments of two degrees) and 
fi ft ft n/ a 
made into a table which la one of the model Inputs. 

' jjjndjt; These are the real and Imaginary parta of the refractive index. As discussed 
earlier la this report, they are used for the determination of Rifl), the Fresnel reflectance. 
Vttuaa for a and k are estimated for the palat surfaces In this study. Moreover, the surface . 
is assumed to be essentially nonconducting so that k * 0. Based on experience with similar 
palat samples, a is taken to be 1.69. For a given sample, a and k can be determined accu¬ 
rately by measuring the Brewster eagle and calculating a and k as outlined in Section 4 on the 
surface model. At the present time the program used, RHOPRIME, does not do this. 

r andO. ^ These parameters are used la the function which provides a correction to the 
program to account for shadowing and obscuration effects resulting from the roughness of the 
surface. Values far r sad A bare been selected, based On observed characteristics of reflec¬ 
tance properties. They hare been established a»r « 15 andQ ■ 40. . 


; Vw. 








7.3. LAMBERTIAN VOLUME MOCEL PARAMETERS 

pj and Py These are the cross components of polarized radiation used in the model to 
account for the diffuse contribution, p ^ - 2p ] , and -2p, a , where p, 5 and/or p,^ is de¬ 
termined by taking the average value of the cross component from the measured data. Ac¬ 
cording to the reciprocity-theorem, p^j - [8®** ®1- ^ is important to remember that p x 

values only used when the volume scatter model Is not used. When the volume model is used. 


7.4. NON-LAMBERTIAN VOLUME MODEL PARAMETERS 

Py This represents the non-Lambertian volume scatter component; it is determined by 
extracting pj ^ or pj at the point which would lie uixier the peak at the zero bistatic scan if the 
measured curve were smooth; p„ = 2p' = 2p‘ , at the peak point. (The fact that a hump some- 

V J, * ' 

times occurs on the measured cross component curve is discussed m the section on Model Vali¬ 
dations.) 




Here again, it Is important to remember that when the Lambertian volume model Is used, 
Py. = 0 and p^ * 0. When the non-Lambertian volume model is used, Py * 0 and = 0. Also, 
Py and are never simultaneously nonzero in models which have been validated to date. 

DP(p), W), Integral parts of SUBROUTINE FUNC, these parameters currently are 

all set equal to 1. They have been Included to provide flexibility far later model modifications. 

7.5. PARAMETERS USED TO GENERATE TYPICAL DATA FOR COMPARISON PURPOSES 
As will be discussed in Section 8, available data representing typical material parameters 
are sometimes useful. For this reason, the model contains a subroutine which can generate 
zero bistatic data, given a suitable set of input parameters. Tae parameters are: o, RPO, Ql, 
and Q2. Normally, however, actual zero btstalic measured data are used, and o = 0; RPO = 0; 
Ql =■ Q2 * 1. ... 
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REFLECTANCE ESTIMATION METHOD 

If one has sufficient information about the target material, reflectance data can be esti¬ 
mated without the use of a computer. In particular, as discussed earlier, the Index of refrac- 
tioo and a fixed bistatic curve are the necessary elements from which to extract the necessary 
parameters. For the sample paints used for the data compilation, the index of refraction is 
assumed to be totally real with n = 1.65 and Y = 0. 

In order to use the above information to generate reflectance data, a set of typical fixed 
bistatic and bidirectional reflectance curves has been generated. These curves are intended 
to simulate a range of paint types from which reflectance values for a particular material 
can be estimated by an interpolation procedure based on fixed bistatic values. Typical curves 
are given in. Appendix IL 

We first find the position af the fixed bistatic for the material of interest relative to that 
of two other fixed bistatics for materials on which we have complete reflectance data. We then 
assume that the same relationship between the three materials will be maintained In the reflec¬ 
tance data. Therefore, by interpolation, we determine the bidirectional reflectance of the ma¬ 
terial of interest relative to the known materials. 

To proceed, we must now choose parameters that characterize the curves and provide a 
basis for interpolation. The parameters selected are: 

(1) the ratio of peak value to that value at which the curve begins to level out or, if it 
begins to rise, the point of minimum reflectance (this is for the like-polarized com¬ 
ponent, i.e., p^or p^). 

(2) the width (In this case, angular distance between peak point and leveling-out or mini¬ 
mum point, of the fixed bistatic curve; for most materials this width is close enough 
to 30° to be as named constant). 

(3) the Lambertian or non-Lambertian character of the material as demonstrated by the 
angular dependence of the cross component in the fixed bistatic scan. 

In the following subsections, we first provtde a step-by-step outline of the interpolation 
proc e d ur e. The tame step-by-step procedure Is then applied In an example worked out tn de¬ 
tail. ■ - 


8.1. PROCEDURE FOR ESTIMATING REFLECTANCE VALUES BY INTERPOLATION 

, (I) Select the actual fixed biataitc curve for the material for which reflectances are re¬ 
quired. 






'Zt Measure ire reflectance at i;>e pea* vatu* cl the curve a,jd the ciicimua v*iue {or 
va-ur at v.r.irn in* carve be?in* to level ou( from (he rvak) Tike the ratio <4 ’hes* tw, values. 

3 • Select two een«ra£ea fixed n'JUiiiC curve* which appear to Bracket i« uicusured curve 
* it tt respect to tr>e ratio described to step (2! 

4i Normalu* ail three carves to tM tiae peas value and determine ute aormaiiiauoe 
factors. 

S) t the measured cure* lies nets era the two reference curves, select an xntie close to 
ai'’ 'he approximate ra.u* for Diaimi or serelirse out f or most materia is and determine the 
puettioa at the measured ixxat between the points at t!>e re a aruftetai curves. 7>># fraction ut 
tne distance between the t»o reference cure* points at which the measured curve pout) lies is 
then takes as aa inierpoia'ioo factor, if This factor will be applied to the reference ret tec - 
tance curves to obtain reflectance values for the material of interest tj the measured curve 
does trot lie between the two artificial corves, the procedure is one of extrapolation, and the 
derived factor will suit be some fraction of the distance between the artificial curve points.) 

■6! for 'hr desired source angle - 2°, 20°. 40°, or 60°), select reference reflectance 
c ur ves from Appendix □ which correspond to the fixed bistatie Curves 

<“7) Beni rdf the reflectance vat see for the r e f e r en ce curves for each receiver f ) ancle. 

f*1 Multiply both reflectance c urves by the asms factor used to eormalue their corre¬ 
sponding fixed bis Ut ICS. 

ft) Using the aster-pots* *ee for extrapnUfMsi factor. If. derived in step (5), find the esti¬ 
mated refhsrisnce rates far each posit toe of v 

Hot# The some renege tor re spo nd to normalised fixed Siststic v .'see. To set the 
abeotd* ref kectsnre vain* we meet divide out the nor waluatine. 

10) 12 nde vs tee* is «ep it tier the ace nsstuing factor • ffeund In step 4) for the material 
rf Inter* *» 

4 1 Aff liCATTCW or fBocrouwt 

We el 11 amuse** fir* 0 f pom* far » . 40*. using the ftxed btstatic for nsmpte A01040 
a# Sh e en in fig. M, Append* L We vtB then enmpnre dm flee esetmeted pomte nttfc the 
scfwl am■ servmeet curve en png* IU at the Dnts Consptlettop|19*4. 2{. The whole p ro cessre 
will he rsrried out ee mdkneO above 

(I) te mpt* A01040 hen bean se l ec t e d and the f u n d 4dm MM curve of fig. 59 will be need 

d, p fpsed) - 0 It 

p (10*) . 9 049 

1.0 It 9 OH - 11« 


90 








<3! Select fixed bistaiic curves for typical niaienals 2 and 3 fuse Ftgs. 63 and 64} with R 
nines that bracket 3.73. ;U u convenient to [.racket the value but cot necessary.) The tvo 
value* can oe used to extrapolate to any value. 

(4 6 St Norraaiiae AO 1640 and material 2 to Use peak at material 3: 


Material 

a 

Measured Value 

NormalLtme Factor 

Marmaiixbd* 


f o a 

0.1* 

5 62 

.3* 

AH16+0 

\30° 

0.035 

5 62 

o.trr 

«* 

f o° 

0 *5 

133 

3.0 

2 

(30° 

0.035 

1.33 

o:*+» 


1 0° 

0* 

1 

.0.* 

3 

130° 

0.38 

1 

0.33 


Far interpolm-on. one hi<hes. itot as r efere nc e: 

p 3 " p is 40 o is - o :rr o tsa , 

ucerpout.o. factor » IF . -7~r- p — * ilJToW * 6331 * 3 554 


Hots that p, refers to p (or material 3: similar swbnenptiac identifies p',^. and p^. 

(61 Far ». » 40°, nee material 2 (Ft*. 79), and material 3 (FI*. 67). 

(7) We ese » r * 20° and 40° is the backer alter had plane, and «■ 0°, 20°, and 40° la the 


forward 


-scatter had plane. Valors shown are (or p (»_, ) where d is 


understood that » * 40 


p(20,0> p(40.0) pf0,0l pfWjaO) 

Material 2: 04)3 0.035 0 04* 0 23 

Material 3 0 32 0 37 0 44 0 *5 


n 


(I) hoc mat i,i at toe I art or (or material 2 wsa 1.33: (or material 3 thta (actor was 1. There- 


pT20,0) p740,0) pT00) pT20,l«» „;•*)») 

i: *t£i ^trr ^.i r 1 


Material 3: 0.32 0 37 0 44 0 05 2.5 

(*) AaMetepfl), pj - < , j**o * * Pj*» p' 1<40 * dj * OTdp^ * Pjltar each pmdtow 

of lat e reel. Therefore: 


) AO1040: 


,*ap •*» sw w 


(10) Dtvidln* the r ene d s «4elepft)by the norsaatUtdf (actor for temple material A01640 
yields the fot t ewM* estimated valeee 



We hare aon obtained Mb tana bidirectional reflectance r a twee at * ( • 40° for Material 
A01*40 of the DWk COMptlMkm where # f « 0* end # f * 20* and 40° ie both the bacfcstaltered 
and l o r w4- «tdiw4d O mtte a t . 


M 





Since measured vaiuea for these reflectances are available in tne Data CcmpUafton [2J, 


.• - w e en s 


*e can nt/ar cocsoare them to our aerived vai'aes to aeterraine how wail the interpolation 
method worts. As page 118 at the Data Compilation snows; 

(.[ '40,0! = 0.06 
p" , (20,0) - 0 042 
p: , (0,0) >0.055 

P ; , ao,i») »o.u 

pj* (40,180) =0.7 

Recall, h ow ev er , that tar trae surface reflectance: 

' V P J 

Therefore, the cross-comaoeeats rear from the same measure m e n t curve are: 

p; # (2O,0) > 0 015 
pf, (40.01 . 0 018 
P t (0,0) -0 017 
pt|{20,180) - 0 018 

p^ <40.180) - 0.02 

Calcuiatlac p * p\ we obtaia *he follflering cocnperisce: 

ii 11 

Messoremewt Estimate 

p, (20.0) . 0 027 0 028 

p( (40.0) • 0 042 0 025 

o', (0,0) -0 038 0.045 

p (20.180) • 0 10s 0 11 

p[(40,180) « 0.68 0 44 

The agreeeiset appears t* be cseeOeet, the larger* discrepancy jaawstutg to sbtart 50%. 


•1 













Appendix I 

FIXED BISTATIC DATA FOR PAINTS FROM 
DATA COMPILATION 

In Section 4.3 on shadowing and obscuration, »* >bowed that Eq. (10) could be used to de- 
rtr* a fixed bUtatlc curve sa a function at source-receiver poeiuon. But by <o doing, oee ob- 
«»*»« exe variation at fixed bistatic curves with source position. Tbs shadowing and abeam 
tide factor described in Sect toe 4.3 xppbes a first correction to tbs basic calculation. 


A fur th er c or re c tion may be applied by averaging values obtained in the first correction. 
Fixed bietatic curves for ?4 rt a mpin s froca the data compilation have been so derived sod are 
included *a this report as Figs. 39 through 43. Each figure includes a listing c t the fixed bi¬ 
static bidi r e c t i o na l reflectances as a function of i ia correspondence with tbs curve. 
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FIGURE «. FIXED BIBTATIC p* FOB 
TYPICAL MATERIAL NO. 3 

FIGURE 8*. FIXED BOTATIC p* FOR 
TYPICAL MATERIAL NO. 4 
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FIGURE 95. p’FOR MATERIAL NO. 4. 
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Appendix HI 

DOCUMENTATION OF BOTRECTSONAL REFLECTANCE PROCSAM (RHOPHIME) 

Program RBOFRIMX is the main calling pr ogr am for sobroatises to read and stare ma¬ 
terials data, perform geometrical calculations, compote bidirectional reflectances far airy 
sauces' receiver position and polarization, and prepare the output in a convenient format. The 
calling sequence, purpose, and calculations performed by each subroutine are given below, fol¬ 
lowed by details on the input data formats. 

DLL DESCRIPTIONS OF SUBROUTINES 

SUBROATTNE DfPATA. This is the first sucroctine called. Material parameters seeded 
for the calculation of bidirectional reflectance are read and stored. Material parameters are 

MAT « material speedier 
* > » « real part at refractive index 

K *>k * imaginary pert of refractive index ■ . ^ 

KX1 mg m diffuse reflectance tar £ poiartxed source 


RE - p x2 - diffuse reflect for I polarized **rc* 
SHOT « p r » votsme reflectance 

liGltA ) . _,i_wt. ._ < . .. _r fm . _ m 


RPO j 

TAO -rfdag) 
OKXGA • u fc 
Qi 
<33 


available to calcalate s' (V« ; Va)** 


EOTC 


i»eg) 

- Q (aegfl 


parameters to calculate a shadowing and obscuration fader to be ap¬ 
plied to p' (cae HHV) is eabrouttae fTTN'C 




«W-« «**) - .„ 

■ A -^. v -: : r ■ ' •< 

ffttl, A title card (optional) Is read \ 


|. table of wt«- 


btstmc Mdir ecdoasl refine* 


to tdautfy oe the printed output tbs cal* 


-t “i* : 

r.. T - 


FACET. The M arti end r eeetver are located ta as earfli-fised. right- handed XTZ ooortH* 
aate syatesa. The XYS eeanpowams of the and normal f a ct o r of tie refltetlag ezrface are read 
(epttoaatL V tbe facet deOaitie* ttrd ts eol suppt i ed. the facet sab normal rector dafanlta to 

t- , •• ; •:<, >/» 

CCMRtTTATBOW WKBTEfT . The epeetftestwa of Marta and re ca tver poettione aad source 
poiarttaflea for coeapwtafiaa of (be fehRrectioeU reflect wee is read. 


• *'Y. *•> 

7:: ' m; v-f...7i- - 7. ' : - w 

:: '.■%? ■ - 7ggIS® 


77-77 






' r ' ; 









2TW = model selector 

TS - zenith ar-gie of source (dee) 

PS = azimuth angle of 3oarce (deg) 

TD = zenith angse cf receiver (dej) 

PD * azimuth angle of receiver ;deg) 

A » intensity o l major axis at polarization ellipse 
B « intensity of minor axis cf polarization ellipse 

PS! = angle at major axis of polarization ellipse from the normal to the plane of incidence 
measured CCW looking into the source, 0 £ PSI s ISO (deg) 

P = polarization d source (0 ; P s 1.0) 
a = handedness d polarization ellipse U1.0 or 0.0) 

MI * material speedier 


SUBROUTINE SCAN . This subroutine defines a sequence d detector positions Icr a speci¬ 
fied source position and polarization. 

CSV « model selector 

TS > zenith angle of source (deg) 

PS « azimuth angle of source (deg; 

TD8 • start zemtn ucgje of receiver (deg) 

TDK - end zenith angle of receiver (deg) 

TSTEP « zenith angle ocas increment (deg) 

PCS « start azimuth angle at receiver (deg) 

PDZ > end azimuth angle cf receiver (deg) 

PS IIP » azimuth angle scan increment (deg) 

A <* I n tens it y eg major axis at polarization ellipse 
B * intensity eg mfeor axle of polarization ellipse 

PSI ■ angle of major axle of polarization ellipse from the normal to the plane d incidence 
measu red CCW lack tag into the source, 0 * PS!« 180 (deg) 

P « polarization of source (Daps 1,0) 
g » h an d s tim e eg poisnxaiioa ellipse (*10 or 0.0) 

Ml - material specifier 

gPOTOUTlNE CIOM. This ssfcroutlae dots the necessary geometrical calculations of 
angle* nee ded tor the bidirectional reflectance calculations (sea Pig. »). 

Oft • 0, 0, 1) U n uaft vector along the earth-fixed Z axis 

PH ■ the angle of the major axis of polarization ellipse from the sormal vector eg the OR, 
t plane measured CCV lock tag into the sourer, 0 * pa s lM (deg) 

_ o»g 

_ . Oft k S 

t *TCKTIT 



















ca v d 

OH x Di 


XA = 




i D * E i • 


XAP = 


YAP = 


OP * E 
[ OP x El 

OP* D 
OP * Di 


C06B = X-D 
C0S3DP = OP-D 
C06BEF = OP • E 
COGBNP = OP • X 

PSIPE = PSI - SIGN(-XAP- OR)AHCOS{XAP • Y) 

= angle of major axis at polarization ellipse from the normal vector of the OP, E 
plane 

PSIDE - PSI - SIGN(Y- D)ARCOS{U • Y) 

* angle of major axis at polarization ellipse from the normal vector cf the D, E 
plane 

WADE > -SIGN(-YA • EiARCCCOXA • YA) 

» angle far transforming the output angle at polarisation from E, 3 plane to OR, D 
plane 

ED PHI > ARCOB(XAP • YAP) • the relative azimuth angle between E and D tn the facet 
coordinate system 

DC = (-6INBEP, 0, COSBZPI » direction of specular ray in the facet coordinate system 
DI * (SINBDP COB ED PH I, SINBDP SIN ED PHI, CQ6BDP) - direction of reflected ray tn 
the facet coordinate system 
NZ1 - DC * OP 
NZ * NZ1 * DC 
DM-DINZ 

PffiBN • 0 IF DM > 0 1 parameter required in TVKC now YUNC 

■ s/2 * ARCOB('DN) IE DM < 0 f for shsdiiwtag tad obecurattan 


SPBKOPTDCB GYRM, GERM does all of the bidirectional reflectance calculations. The 
sartor oat las requires: 

F • series of switches which can be set \r . A on \l) to reduce the anmtoer of redundant compu- 
tafioas whee GYRM Id used as part ft > multifaceted target model 
COBS - defined la SUBROUTINE GEOM 
COBBDP - defined tot 3UBRC* INS GXOM 
COBBZP - defined la SUBROUTINE GEOM 



'iwmm 


- - -r- ' ? • is/y fiW- ;.:r. v ■ 
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1 


j 


j 

i 


CCSBNP = defined in SUBROUTINE GEOM 
P3IPE = defined in SUBROUTINE GEOM 
PSIDE = defined in SUBROUTINE GEOM 
WADE = defined in SUBROUTINE GEOM 

AP = area d facet (if AP * zero, GERM returns a bidirectional reflectance Stokes rector; 
if AP * 0, GERM returns a Stokes vector for the reflected radiant intensity for unity 
irradiance in the incident beam) 

MX - material specifier (available in COMMON) 

ISW x model selector (available in COMMON) 

W = Ta re length specifier (available in COMMON), not used 

TABUS = array containing all of the materials properties data read in SUBROUTINE IN- 
DATA 

GERM returns the bidirectional reflectance Stokes vector (AP = 0) or radiant intensity Stokes 
vector (AP #0). 


Ill * Stakes vector for surface pins Lambertian model with polarized scarce 
B1 * Stokes vector for surface pins Lambertian model with unpolarized source 
US * Stokes vector for non-Lambertian volume model with polarized source 
B3 * Stakes vector far nan-Lambertian volume model with unpolarized source 
114 « Stokes vector for combined model with polarized source 
B* * Stakes vector for combined model with unpotawtzsd source 

FUNCTION GETDAT returns the appropriate material parameters for bidirectional re¬ 
flectance calculations, namely N, K, RX1, KX2, RHOV, RCOSBNP, DPO, DP90, F, a 

FUNCTION FUNC prrbrides the optional capability for deriving RCOSBK? analytically (if 
SIGMA * 0)and far deriving a shadowing and cbecaratlcn correction factor (optional) to the 
RCOSBNP used in the specular model. In addition, the depolarization factors DPO(B) and 
DP90(B), an well aa F(B) and C*BNP) nee de d in the volume model, are defined analytically, 
FUNCTION FUNC currently yields 

DP0(B) • U> 

. . DF*KB}» 1.0 
F(B) <*1.0 



* / 



i 








The sn&doung and obscuration {actor applied to RCOEBNP (measured values read during the 
input phase of RKOPP.IME or detir.ea analytically in FITNC) is 


i BNP TAU 
* OMEGA e 


OMEGA 


t PHIEN _ BEP 
OM£ G A OME G A 


(a) suhfacs-plus-ia.mbertun model c alculation 


N + 1 r - K 2 [V2 - CCSB) 2 - V3 

KU - o 1 '} 


(N - 1)“ - K" (V2 1- COSB)“ * V3 
normalized reflectance for 1 poiarized incidence 


(V2COSB + COE”B - l)“ * V3CCS^3_ 

rww - --^^^ 

(V2COSB - COST'S + 1} + V3COS“B 
= normalized reflectance fcr I polarized incidence 


V2 * 'H 2 - K 2 - V > CCS 2 Bl 2 ■ (N 2 - K 2 - 1 + COS 2 B) 

vn l^^K 2 * (N 2 - K 2 - 1 ♦ COS 2 B) 2 - IN 2 - K 2 - 1 COB 2 B) 

V3 * 2 

C B - 0 the calculation is made for a plane polarized source (polarization angle PSI). The 
calculation Ignores the induced elliptical polarization for K < 0. 


PEED =. A TAN 


TAN PBXDE • SION (COSATAN(N) - COBB) 


> polarization angle with reaped to D, E reference plane, after reflection 

C > 1 if AP * 0, then a bidirectional reflectance Stakes vector la computed 
C «- AP • COSBEP * COSBDP If AP # 0, than a reflected radiant Intensity Stokes vector Is 
competed 

111(1) - C SillirKiffinp <*° • COsS'SnjK ♦ R80'SD^PfflDE 

111(2)-C ^ ^p^ fi e BDp (*» • C03 2 PSIDE tRWSIN 2 PgffiEJCOSa (PSIED-TVADE) | 

• J 

111(3) - C cO^bIp^& BDP * COg2 PS1DI ♦ PSIDE)SIN3 (PSlED-WADE)j 








m(1 > ■ C [j §BI%C65i)P l (R0 - RS0 > 


J {tan. + RX2) 


Ei(2) = c h m - E90) c< -^ ( - v/AUE) 

m m 

m(3) = C iffiiSlc6BDp| (RC - «»> 5012 «' WAI3E ' 


g H * *1 the calculation Includes the phase difference and elliptic ity induced by reflection 
for K * 0 and is an exact treatment of the Fresnel equations. 

SUBROUTINE ELEP3T (AA, AB, PSIDE, H;.AA1, AA2, D) defines the following input ellipti¬ 
cal polarization parameters: the amplitudes perpendicular (AA1) and parallel (AA2) to the D, 

E plane and the relative phase D = <p» - $1 of the amplitudes of the major (AA) and minor (AB) 
axes of polarization ellipse; the orientation of the ellipse with respect to the D, E plane; 

PSIDE; and the handedness, H. 

DR * 4t - *i induced by reflection 

FOR K = 0, DR * 0 if COSB < COS ARTAN(N) 
if COBB > COS ARTANCNJ 


FORK * 0, DR ~ s * ATAN 


( 2 - CQS 2 B)C06B 

l (1 - COsS} 2 - COS 2 B(V 


tf ( )<0 


DR- - ATAN 


I 2 VvTa - COB*B)CCBB -__V 
^ (1 - CCB^BT - COe^B(V?~*- V5)y 


If n > a ■' 

The intensities AIR and A2H and the relative pease of the parallel and perpendicular compo¬ 
nents cf the reflected radiance induced by the reflections, are 

A1S >A1'BD .. '£.? " . . 

A2R-A2-S0 '■ ' ~ _ . y* ' 

dr-dr+ D ■ : > : r '7.:X4y 

SUBROUTINE EUPBZ (AA1, AA2, DR, AAR. ABR, PSCED, HR) defines tbs eUlpticaUy 
pitarized r«fk..-ted radiance as amplitudes AAR and ABR of the major and min o r faces, tbs 
ssfte at the ellipse relative to the D, E (dans, PfflXD, and the handedoMs, HR. PfRDB • 
PSZKD-WADE is tbs angle that tbs polarization ellipse of the reflected radiance mafcss wttb 
tbs norms! veetor tc (be OR, Dplan*. • • . ’ v 

CHI - HR' ATANfABR/AAJl) is the parameter used to define tbs clltptictty at the reflected 


radiance. 

C-ltf AP-0 


• • i?-r-- .... - 










i 


! 

i 


V- 


Si 

w 


U 

f, 


C = AP ■ CCSBEP • COSBDP if AP * 0 
111(1) = C| 

111(2) = Cj 


111(4) = C 


FAR + BR 

RCOSBNP 

L A +■ B 

COSBEP COSBDP 

Far ^ br 

RCOSBNP 

L A+ B 

COSBEP COSBDP 

Far - er 

RCOSBNP 

L A +■ B 

COSBEP COSBDP 

AR +■ BR 

RCOSBNP 

L A +• B 

COSBEP COSBDP 


- (RX1 COS 2 
COS(2PSIDE)COS(2CHI) 


(PSIPE) + RX2 SIN 2 (PSIPE)] 


SIN(2CHI)| 


121(1) = Cl 


RCOSBNP 


COSBEP COSBDP 2 


1 (RO + R90) * |(RX1 + RX2)| 


121(2) = C 


| RCOSBNP 


[COSBEP COSBDP 2 


i(R0 - R90) COS(-2WADE) 


121(3) = C 


; RCOSBNP 


4 (RO - R90) SIN(-2WADE)| 


mCOSBEP COSBDP 2 

121(4) = 0 

(b) VOLUME MODEL CALCULATION 

The angular-dependent, volume reflectance model, Stokes vector is given by 
I13(1) = C bP90 7 l>DPO) ISpSjSbD pt 00 ® 2 PSIDE-DP90(UDP0) * 
SIN 2 PSIDE • DP90(l+DP90)j 


113(2) = C 


Sp^ R npl COs2 PSIDE-DP90O-DP0) 


DP90( 1+DPO) C06BEP+C06BDP 
SIJ^PSIDE • DP0(l-DP90)]COS 2 2AD 


113(3) = C 


2RHOV-FG 


[COS 2 PSIDE ■ DP90(1-DP90) 


DP90(l+DP0)CO6BEP+COSBDP 
Stt^PSIDE-DPOd-DFWjSIN 2 2AD 

I23 * 1 * “ C DP90(1+DP0) COSBEP+COSBDP 2[ DP9 °( 1+DP °) + DPO(1+DP90)] 
= C DP90(1+DPO) COSBEP+COSBDP 2 f DP90 * DP0 ) COS(-2WADE) 
123(3) = C Dp9o(i + 0po) COSBEP+COSBDP 2 f DP9 ° “ DP0 I SIN(-2WADE) 


where 

C = 1 for AP = 0 

C * AP • C06BEP • COSBDP for AP * 0 
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The angle of polarization d the reflected radiance, AED, from the normal vector of the D, E 
plane is 


AED = ATAN 


/ PP0(1 - DP90) 
DP90(1 - DPO) 


TAN(PSIDE) ■ SIGNtCOS ATAN(N) - COBB) 


and the angle at polarization referred to the OR, D plane is 


AD = AED - WADE 


SUBROUTINE OUTPUT . This subroutine prints the Stokes vectors for the bidirectional 
reflectance (AP = 0) or reflected radiant intensity for unit Incident irradiance (AP * 0) for the 
surface model, the volume models, and for the combined specular and volume model. Stokes 
vectors are printed far a completely polarized source,for a completely unpolarized beam, and 
also for a partially polarized beam (polarization defined by the input parameter P). 

In addition, several calculations are made with the Stakes vectors. For a bidirectional 
reflectance (or radiant intensity) Stakes vector, the bidirectional reflectance (or radiant in¬ 
tensity) for a receiver polarized 1 or 1 to the OR, D plane is 

A + n 

receiver ^ j— 


receiver j = —^— 

where the Stakes vector is at the form: 


The angle at the major axis of the reflected radiance and the percent polarization at the 
reflected radiance are also given; they are 

. j 1 tm.w iC| on < AT < nr, (looking into the source, AL > 0 Is a CCW angle; 

AL *2 ATAN !bI - 90s AL-90 al < o u E cw angle) 



The output Includes TS, PS, TD, PD, P, as well as the input and output values of A, B, PSI, H 
from the surface model calculation (If the Input H ■ 0, the Input and output values of A, B, H 
default to 1, 0, and C). 

SUBROUTINE ELIPS1 (A, B, PST, H; Al, A2, DELTA) . The basic equations which relate 
two specifications of an eillpttcally polarized beam (A, B, PSI, H) and (Al, Bl, DELTA) are 

tan a * Al/A 2 0 s a * v/2 

tan x * *8/A for * - »/4 * X * »/♦ _ 
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tan 2v = tan 2 a cos 6 
sin 2\ = sin 2a sin fl 


from which we obtain 


7 2 

2 - sin 2x * tan 2’*' 

sin 2a = - 2 - 

1 ♦ tan 2v 


or equivalently 

cos 2a = cos 2x cos 2^ 

Subroutine ELIPS1 determines Al, A2, DELTA from A, B, PS I, H 


LAMBDA = \ A 2 + B 2 


If B - 0, A1 - LAMBDA COS PSI 
A2 = LAMBDA SIN PSI 
DELTA - 0 when 0 s PSI s # 2 
DELTA - t when tr/2 * PSI * t 


Otherwise: 

CHI = H • ATAN(B/A) 

n = COS 2CHI COS 2 PSI I 

ALPHA = 1/2 ARCOS(-TI) II t/4 s PSI s 3r/4 

* 1/2 ARCOS(TI) U PSI < »/ 4 or > 3#/4 
II ALPHA = 0, A1 = LAMBDA A2 = 0, DELTA = 0 
II ALPHA * */4, Al = A2 = LAMBDA/V?, DELTA = 2CHI H PSI = »/4, 

= H-r - 2 CHI U PSI = 3ff/4, 

D ALPHA = s/2, Al = 0, A2 « LAMBDA DELTA = 0 
Otherwise 

TI = ISIN2CHI/SIN2ALPHA! 

MU * ARSIN TI 
Al » LAMBDA COB ALPHA 
A2 » LAMBDA SIN ALPHA 
C06D = TAN 2PSI/TAN 2ALPHA 
Ef COBD > 0 DELTA » H MU 
H COBD < 0 DELTA - H-(* * MU) 

SUBROUTINE ELIPS2 (Al. A2, DELTA A B, PSL Hj . Subroutine ELTPS2 determines A 
B, PSI, H Irom Al, A2, DELTA 

M 





I 

i 


LAMBDA = VaI 2 + A2 2 

If Al = 0 or A2 = 0, then A = LAMBDA, B = 0, H = 1, and 
PSI = 0 if A2 = 0 
PSI = x/2 if Al = 0 

a A1 = A2, then CHI = 1/2 IDELTAi, A = LAMBDA COS CHI, B - LAMBDA SIN CHI, and 
H = 1 tf DELTA > 0 
= -1 tf DELTA < 0 
PSI = »/4 if CHI < J/4 
= 3x/4 a CHI > 5/4 

H DELTA = ± 5 , A = LAMBDA, B = 0, H = 1, PSI = » - ATAN ^ 

H DELTA = 0, A = LAMBDA, B = 0, H = 1, PSI = ATAN 

a DELTA = ±x/2, H = +1 a DELTA > 0 
-l a DELTA < 0 
H A1 > A2, A = Al, B = A2, PSI = 0 
B A1 < A2, A = A2 r B * Al, PSI - »/2 

Otherwise 

a Al > A2, ALPHA = ATAN A2/A1 

CHI * 1/2 ARSINiSIN2ALPHA SINDELTAl 
LAMBDA » i TAN2ALPHA C06DELTAI 
A - LAMBDA COS CHI 
B - LAMBDA SIN CHI 
H - 4 a DELTA * 0 

Part I; 0 < I DELTA I < »/2; PST * 1/2 ATANLAMDA 
Part 2: x/2 < I DELTA! < x; PSI « x - 1/2 ATANLAMDA 

and 

ff Al < A2, 0 < iDELTAi < »/2 PSI - x/2 - 1/2 ATANLAMDA 
x/2 < IDELTAi < x PSI * x/2 + 1/2 ATANLAMDA 

HL2. INPUT DATA FORMATS 

The inpat to the RHOPRIME program la segmented Into logical blocks. Each block is Ini¬ 
tiated by a block header and terminated by an end card. Blocks may be in any order, but a data 
block is assumed to precede any computation request blocks or scan request blocks. If a block 
header specifies an invalid block types, all input up to and including the next end card is Ignored 

DATA TABLES BLOCK. The data tables block specifies all physical characteristics of the 
materials to be studied. The block header ts one card with the following format: 


i 



/ 








Columns 


Descrintion 









T^r-^^-^r. *&&&**?**: sw 




Columns 


Description 


’COMP’ 
ignored 
model selector 
ignored 


The model selector, ISW, is: 

1—if specular and diffuse models are desired 
3—if volume model is desired 
7 —if combined model is desired 


Following the block header, computation requests are processed sequentially until an end 
card is encountered. The format cf a computation request is: 

Columns Description 

1-4 blank 

5-9 ignored 

10-16 source zenith (deg) 

17 ignored 

18-24 source azimuth (deg) 

25 ignored 

26-32 detector zenith (deg) 

33 Ignored 

34-40 detector azimuth (deg) 

41 ignored 

42-48 polarization major axis length 

49 ignored 

50-56 polarization minor axis length 

57 Ignored 

58-64 angle of source polarization (deg) 

65 ignored 

68-72 source percent polarization + 100 

73 Ignored 

74-76 handedness at polarization (if * 0, elliptical polarization is assumed) 

77 ignored 

78-80 material index 

SCAN REQUEST BLOCK . If a scan at the detector zenith and/or azimuth is desired, a 
scan request block may be used. The block header Is one card with the following format: 


Cotnmns Description 

1-4 SCAN' 

5-19 Ignored 

20-25 model selector 

26-80 ignored 

One card follow* the block header giving all required parameters. The format at this card 
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■ VT . y it 




wrap* 


••*■’ 4.1*. i 


Columns Description 

1-6 
7-12 
13-18 
19-24 
25-30 
31-36 
37-42 
43-48 
49-54 
55-60 
61-66 
67-72 
73-76 
77-80 

TITLE SPECIFICATION BLOCK . A title may oe printed at the top of each page of long 
form output using the title specification block. The block header is one card in the following 
formal: 

Columns 


1-4 
5-19 
20-25 
26-80 

One card following the block header specifies the title. The format of this card is: 

Columns Description 

T-60 title 

61-80 ignored 

FACET DEFINITION BLOCK . If default facet definition is not desired, the facet may be 

redefined using the facet definition block. The block header is one card In the following for¬ 
mat: 


Description 

'TITL' 

ignored 

blank 

ignored 


source zenith (deg) 
source azimuth (deg) 
initial detector zenith (deg) 
final detector zenith (deg) 
zenith increment (deg) 
initial detector azimutn (deg) 
final detector azimuth (deg) 
azimuth increment (deg) 
polarization major axis length 
polarization minor axis length 
angle of source polarization (deg) 
source percent polarization + 100 
handedness of polarization 
material index 


Columns 


Description 


1-4 

5-19 

20-25 

26-80 


'FACE' 

ignored 

blank 

ignored 


One card following the block header defines the facet. The format of this card Is: 


Columns 

Description 

1-4 

blank 

5-9 

ignored 

10-16 

facet area (default « 0) 

17 

ignored 

18-24 

facet normal - x (default ■ 0) 

25 

ignored 

26-32 

facet normal • y (default * 0) 

33 

Ignored 

34-40 

facet normal - t (default * 1) 

41-80 

blank 


• ( 

i 
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Appendix IV 

INSTRUCTIONS FOR USE OF FROCRAM 
WITH SAMPLE COMPUTER OUTPUT 

The program documentation in Appendix HI, together with the sample computations in¬ 
cluded in this appendix should enable the user to (1 > modify this program to accommodate the 
requirements of his own computer and (2) verify output from his modified program by com¬ 
parison with the samples given herein. 

Note that the input parameter values shown in Table HI are the ones with which the program 
has been run. 

Sample outputs presented in this appendix include 

(1) a listing of the input information (Table III) 

(2) the computed output of the program (long form) (Table iv) 

(3) a short form of the computed output, containing only that information necessary to feed 
into a computer program for the purpose of obtaining plots of the data (Table V) 

The three tables mentioned above appear at the end of this appendix. All of the sample in¬ 
formation is keyed and labelled so that elements may be identified easily. However, the further 
descriptive detail below may be helpful in studying the samples given. 


KHOPRIMS Input Listing 

The following items appear across the topof Table IH. Or. line 2: 

o • real part of index of refraction 
k » imaginary part of index of refraction 
p xl * cross component <2^ )1 ^ {Qr gurlace model 
- cross component (2pjj ) 
py « volume component used for volume model 
SIGMA ■ generating function parameter 
RPO x generating function parameter 

And on line 3: 

r x chad owing and obscuration parameter 
0 x shadow lag and obscuration parameter 
Qt • generating function parameter 
Q2 ■ generating function parameter 


tabulation which, in thla 


following these Items tn Table in Is the p'fc # ; 8 .* \cos S 6 tabulation which, in this 

Van fi It/ 2 

case, was extracted from measured data and determined from the xero blstatic scan. Alterna¬ 
tively, such s tabulation can be generated by use of a generating function specified in the SUBROU- 


ttnx ruse. 








Note ia the sample input information ot Table HI that vahisa are provided for p^, and p^ n 
and also for p^. In practice, p^, and p^ ‘**11 be used or py will be used; ail three values 
will never be nonzero simultaneously. 

If the table is supplied as part of the input, the parameters SIGMA and RPO are set to 
0 and Q1 * QZ * 1. 

The p' /ft S \c°* 2 8 . tabulation is followed by scan request information telling 

I A A A A I A 

Von n a/ a 

the computer what source-receiver combinations are to be computed and what model is to be 

selected; 

ft * ft for source 

♦j * 4 for source 

8 f ^ * Initial ft for receiver 

9 • maximum ft for receiver 

r2 

® r j « size of angular step for ft^ scan 
♦ r l = # for receiver 

« e for receiver (value for second scan) 

♦r3 » size of angular step for * r 
A * semi-major axle of polarization ellipse (normalized to 1.0) 

B * semi-minor axis of polarization e.lipna (B * 0 implies linear polarization) 

PSt » softs at source polarization 
P • per cent polarization (1.0 » 100%) 

MI * material index 

BV * T for combined model. (When volume modul Is used, set j » p^ * 0.) 

Note that in addition to these input parameters, others must bn sided in the SUBROUTINE 
FUNC: 

DF0, DP90 » depolarUations for p er pen di cular sad parallel c ompo n ents at inc iden t beam 
f, g ■ votame modal parameters. 

Far dm sontertals in dm sampln listing, vaknee for DP0, DFtO, £, sad % ken been set equal 

to L0. t ;• • •' . . • A*-..- v.’.~VV:..,. 


As exem pt tfl ed by Table IV, each paga of dm computed output corresponds to oue source- 
receiver configuration. Bams at dm upper left arc self-explanatory. H owe v er tt should bn 
borne in mind that MAJOR refers to dm semi-major elliptical axis (a), which is take* to be 
Ml. fines MINOR, which refers to dm semi-minor axis (b), is 0, the MAJOR-MINOR cos*feta¬ 
tion implies Baser polarisation with polarisation angle Pfi far the incident bnun HANDED • 

0 whenever the polarisation is liaaar only. 


: * * r *■ ■' ykiatjSr • 

• <- ?.. - f-t- l - • V.’v v.‘ v ;.•» 


, Mb » 










The entries in the three main columns are reflectances From the top, the first four 
entries in each column are the surface niotlel elements of the Stokes vector which describes 


the polarization state of the beam as it leaves the target: 

A * total reflectance 

B * reflectance «vith receiver polarization angle * 0 (perpendicular polarization) 

C - reflectance with receiver polarization angle = 45° 

D » reflectance with receiver circularly polarized 

The second four entries, still In the surface model block, are 

4^.0 

—~ - reflectance recorded from receiver with analyzer set for perpendicular polarization 

A-n 

— y ~ - reflectance recorded from receiver with analyzer set for parallel polarization 

AL - angle of polarization for reflected beam 
P * percent polarisation of reflected beam 

Thus far the first two blocks of four entries have been discussed. The foregoing, as pre¬ 
viously stated, apply to the surface plus Lambertian volume model. 

The third and fourth blocks apply to the non-Lambert tan volume model and are to be inter¬ 
preted In exactly the same manner as above. 

The fifth sad sixth blocks consist of the sum of the surface + volume models and are printed 
out ten convenience. 

Note that in the volume model output and la the summed output, item D (ctrcularly- 
pcls. tzed component) is not present. 

Computer Output (Short Form) 

The short form of the computer output consists of the information in the last four entries 

A-B 

of the summed otkput (surface * volume), -j-, , AU P (see Table V). Moreover, the 

data are compressed so that, whereas the long form has only one source-receiver configuration 
per page, the short form cots sine a complete seas is one block. 

One seas consists of four Item numbers. Preceding each of the first two item numbers in 

each seta are 

Wavelength (1.0* pm) 

^( 0 °) 

♦ r <Q°or 190°) 

1 

The first item number in each scan coots las —y—. Each output entry is precoded by the 
teas angle—i.e^ 0.0, 0.0290 moans that me reflectance at 6 * 0 is 0.0268. The second Uem 
number contains The third item number contains the polarization angle, AL, at each 

receiver angle. The fourth item number can tains the percent polarization. 

The scans are is the tarns overall order as those in the long form of the output 

*$■ 
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Appendix V 

RHOPRIME PROGRAM LISTING 
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» M np»i**E *s op nj.in.Tj 

1 0T*ER3XnN *(590),npf3),P(1).O(ii,n8(3),t.APELC15i).TA*Lr(E0n) 

2 _ _ £ QUIV *Lrnr£ (TA4 lF,«) ___ 

J ~ irtegpr cnop,T**t.ES/ , T*m. , /»cnHP/'cnMP'/,ERn/'ERo •/,scam./'scar*/“ 

« IMTFGFR TTTlE/'TITL’/.FAtET/'PAfE'/ 

5 RPAL T2l(4).111(41.T23<3).113(1),124(3),114(1) 

b C n H"OR «I.13*,*,TABLE,I?l.Ill,123,113,I?«,X1« 

7 COHMOM /C>»PT/P31,P01,4ETA.8fTAB.RPa.CnslNe,3I6HA,PHlEM,4EP,T3 

_8_0*T» ( IP/Q. O .a.O.l.Q/.AP/O.o/.nR/O. Q,Q. 0,1.0/ __ 

9 C»***•••«**•»••••* 

10 c 

tl C FOR“AT3 


FnRRAt(A«,lSX.I6)__ 

FnB*AT(A4,S*,*(E7.2»1»),P3,e,iX,13) 

pfl(J*»T( • J909MAI. TPRMINATION * ) 

rCR* aTC• 1*, 1 •***• ERO-OP-PILE ERCrtURTPRED* ) 

pnflMATClTAbLE REAd PR»C» — COn^ITIDn CODE -*.F3. 

FHbmaT( • 1 •, • *.«•• harming £np IR COMPUTATION REQUESTS') 
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